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Abstract:  Verwey  has  suggested  that  the  low  temperature  transition  in  mag- 
netite is  due  to  the  ordering  of  the  ferrous  and  ferric  ions  in  the  octa- 
hedral interstices  of  the  spinel  lattice  and  has  proposed  a specific 
ordered  structure.  This  arrangement  would  require  that  the  crystal 
symmetry  changes  from  cubic  to  orthorhombic.  X-ray  diffraction, 
electric  conductivity  and  magnetization  measurements  confirm  that  the 
transition  leads  to  an  orthorhombic  structure.  The  c axis  of  this  ortho- 
rhombic structure  is  established  along  the  cube  edge  closest  to  the  di- 
rection of  magnetization  in  each  domain.  An  external  magnetic  field 
applied  to  a crystal  of  magnetite,  as  it  is  cooled  through  the  transition, 
establishes  a preferred  orientation  for  the  c axis  throughout  the  whole 
crystal  by  aligning  the  domains  akove  the  transition.  The  application 
of  a sufficiently  strong  magnetic  field  below  the  transition  can  cause  this 
c axis  to  switch  to  a new  direction.  This  process  involves  a co-operative 
rearrangement  of  the  ferrous  ions  in  new  sites,  and  is  accompanied  by 
relatively  large  changes  in  dimensions.  In  stoichiometric,  synthetic, 
single  crystals  the  transition  occurs  at  119. 4°K  and  is  marked  by  an  ab- 
rupt decrease  in  the  conductivity  by  a factor  of  90  in  a temperature  in- 
terval of  1°.  No  thermal  hysteresis  is  observed.  The  conductivity  of  a 
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crystal  cooled  through  the  transition  in  a strong  magnetic  field  is  aniso- 
tropic below  the  transition.  It  can  be  represented  by  the  equation 
<r  = A + B(1  + cos 2 8),  where  8 is  the  angle  between  the  c axis  and  the  di- 
rection of  measurement.  The  ratio  B/(A  + B)  increases  rapidly  as  the 
crystal  is  cooled  from  the  transition  to  90°K,  indicating  a progressive 
increase  in  the  long-range  order.  The  magnetic  anisotropy  energy  be- 
low the  transition  has  the  form  consistent  with  orthorhombic  holohedral 
symmetry.  The  c axis  is  the  direction  of  easy  magnetization  below  the 
transition,  and  the  anisotropy  energy  is  very  much  larger  below  than 
above.  The  anisotropy  constants  have  been  determined  at  85°K  from 
measurements  of  the  areas  between  the  magnetization  curves  in  differ- 
ent crystallographic  orientations.  If  the  orthorhombic  axes  are  differ- 
ently oriented  indifferent  regions  of  the  crystal,  the  magnetic  be- 
havior of  magnetite  becomes  still  more  complicated. 


Magnetite  has  been  intensively  studied  in  this  laboratory  in  recent  years, 
both  as  a prototype  of  the  ferrites  and  because  of  its  unusual  behavior  at  low 
temperatures  (from  -140°  to  -190°C).1’  ^ The  initial  permeability  passes  through 
a maximum  at  ca.  -140°C  and  then  decreases  rapidly  to  a very  small  value  be- 
low -155°C.  Bickford*)  has  shown  that  this  maximum  is  due  to  the  vanishing  of 


the  crystalline  anisotropy  energy.  At  ca.  -155°C  a transition  is  observed.  The 

3) 

specific  heat  traverses  a maximum: 4 * *  7 the  electrical  conductivity  decreases  by 

4) 

a factor  of  ca.  100;  and  the  material  is  much  harder  to  magnetize  to  saturation 


1)  L.  R.  Bickford,  Jr.,  Te-h.  Rep.  23;  Phys.  Rev.  78,  449  (1950). 

2)  C.  A.  Domenicali,  Tech.  Rep.  25;  Phys.  Rev.  78,  458  (1950). 

3)  R.  W.  Millar,  J.  Am.  Chem.  Soc.  51,  215  (1929). 

4)  E.  J.  W.  Verwey,  Nature  144  , 327  (1939);  E.  J.  W.  Verwey  and  P.  W.  Haayman, 
Physica  3,  979  (1941). 
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below  the  transition.  7 The  direction  of  easy  magnetization  below  the  trans- 
ition is  influenced  by  a magnetic  field  applied  as  the  sample  is  cooled  through 

the  transition. 

4) 

Verwey  7 proposed  that  above  the  transition  magnetite  has  the  inverse 
spinel  structure,  and  that  the  transition  itself  is  due  to  an  ordering  of  the  ran- 
dom distribution  of  ferrous  and  ferric  ions  which  occupy  the  octahedral  inter- 
stices in  the  crystal.  The  ordered  arrangement  shown  in  Fig.  1 appeared  to 

7) 

him  as  the  most  likely  one.  7 Verwey  suggested  that  the  c axis  of  this  ordered 
structure  would  lie  along  the  cube  edge  closest  to  the  direction  of  the  external 
field  applied  during  cooling  through  the  transition. 

The  object  of  the  present  research  was  to  study  this  unusual  type  of 
transition  by  X-ray  analysis  as  well  as  by  electric  and  magnetic  means.  This 
investigation  was  undertaken  on  the  single  crystals  of  magnetite  grown  by 
Smiltens®)  in  this  laboratory. 

CRYSTAL  STRUCTURE 

The  ferrites  have  the  same  crystal  structure  as  spinel  (MgOAl^O^).  The 
unit  cell  of  the  structure  contains  32  oxygen  ions  in  an  almost  cubic  close-packed 
arrangement.  The  metallic  ions  occupy  two  types  of  interstices  in  this  oxygen 
lattice:  8 tetrahedral  and  16  octahedral.  In  simple  ferrites  two  variations  can 

t 

occur.  In  those  with  a normal  spinel  structure,  one  kind  of  ion  occupies  each 
type  of  site.  In  the  inverse  spinel  structure,  the  octahedral  interstices  are 
occupied  by  two  kinds  of  cations,  distributed  at  random  over  the  16  sites. 

5)  P.  Weiss  and  R.  Forrer,  Ann.  phys.  £lo]J  12 , 279  (1929). 

6)  C.  H.  Li,  Phys.  Rev.  40_,  1002  (1932). 

7)  E.  J.  Verwey,  P.  W.  Haayman  and  F . C.  Romeijn,  J.  Chem.  Phys. 

_15,  181  (1947). 

8)  J.  Smiltens,  Tech.  Rep.  49,*  J.  Chem.  Phys.  20,  990  (1952). 
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Neutron  diffraction  measurements  have  established  that  magnetite  has  the  in- 

9) 

verse  structure. 

An  X-ray  investigation,  the  details  of  which  have  been  published  else- 
where, showed  that  the  symmetry  of  magnetite  below  the  transition  is 

orthorhombic,  with  unit  cell  dimensions:  a = 5.912,  b = 5.945,  c = 8.  388A 
at  78°K.  The  c axis  lies  along  one  of  the  original  cube  edges,  and  the  a and 
b axes  lie  along  the  face  diagonals  normal  to  the  c axis.  These  results  are 
in  agreement  with  those  obtained  by  Bickford^  from  measurements  of  the  de- 
formations of  single  crystals  on  cooling  through  the  transition.  This  is  also 
the  unit  cell  that  would  be  expected  on  the  basis  of  the  ordered  arrangement 
of  Fig.  1. 


ESTABLISHMENT  OF  THE  MAGNETIC  AXIS 
The  orthorhombic  axes  can  have  any  of  six  possible  orientations  relative 
to  the  cubic  axes.  In  the  absence  of  any  external  conditions  to  make  one  or 
more  of  these  orientations  preferred,  they  should  all  occur  in  approximately 
equal  amounts.  The  effect  of  a magnetic  field  in  making  one  or  more  of  these 
orientations  preferred  was  studied  using  the  null -coil  pendulum  magnetometer 
developed  in  this  laboratory  by  Domenicali.  Some  modifications  were  made 
in  the  magnetometer  to  permit  its  use  down  to  liquid  helium  temperatures.  The 
case  of  the  instrument  was  evacuated  and  helium  gas,  admitted  through  a needle 
valve,  was  used  to  regulate  the  rate  of  warming  by  controlling  the  heat  loss  from 
the  sample.  The  sample  was  mounted  in  an  aluminum  boat  in  good  thermal  con- 
tact with  the  null  coil,  which  also  served  as  a heat  source.  The  measurments 
were  made  on  a circular  disk  cut  from  a (lOO)-plane  slice  of  a single  crystal  of 

9)  C.  G.  Shull,  E.  O.  Wollan  and  W.  C.  Koehler,  Phys.  Rev.  84_,  912  (1951). 

10)  S.  C.  Abrahams  and  B.  A.  Calhoun,  Tech.  Rep.  61. 

11)  L.  R.  Bickford,  Jr.,  Revs.  Mod.  Phys.  25,  75  (1953). 

12)  C.  A.  Domenicali,  Tech.  Rep.  24;  Rev.  Sci.  Inst.  21,  327  (1950). 
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magnetite.  Tne  edges  of  this  disk  were  ground  by  hand  until  it  was  approxi- 
mately an  oblate  spheroid  with  axes  of  4.  55  and  0.  89  mm. 

The  smallest  magnetic  field  that  sufficed  to  establish  the  magnetic  axis 
was  determined  by  cooling  the  (100)  sample  through  the  transition  with  mag- 
netic fields  of  various  strengths  applied  parallel  to  the  Of)  direction.  The 
sample  was  then  demagnetized  and  the  magnetization  curve  in  the  [OOI]  di- 
rection measured  at  -195°C.  From  the  curves  shown  in  Fig.  2 it  was  clear 
that  the  minimum  field  necessary  to  completely  establish  the  magnetic  axis  is 
ca.  1000  oersteds,  and  that  this  is  approximately  the  same  field  needed  to  pro- 
duce technical  saturation.  The  magnetic  axis  is  established  along  the  cube  edge 
nearest  the  direction  of  magnetization  in  each  domain.  The  function  of  the  ex- 
ternal field  is  to  align  the  individual  domains  so  that  throughout  the  sample  there 
will  be  a unique  magnetic  axis.  It  is  obvious  that  this  magnetic  axis  must  be  the 
c axis  of  the  orthorhombic  structure. 

When  a sample,  which  has  been  cooled  in  a field  to  establish  the  c axis 
throughout  in  the  same  direction,  is  demagnetized  the  c axis  is  the  direction  of 
easy  magnetization  and  the  domain  pattern  in  the  demagnetized  state  will  consist 
of  antiparallel  domains  aligned  along  the  c axis.  This  was  verified  by  measuring 
the  longitudinal  magnetostriction  (see  Appendix  ?. ) in  the  £lOo[]  bar  after  cooling 
through  the  transition  in  a field  parallel  to  its  length  and  then  demagnetizing  the 
bar.  The  longitudinal  magnetostriction  in  fields  up  to  ca.  10,  000  oersteds  was 
extremely  small  (6 l/t<  10  ) indicating  that  the  magnetization  in  this  direction 

was  almost  entirely  due  to  the  motion  of  180°  walls. 

In  additional  experiments  the  (100)  sample  was  cooled  through  the  trans- 
ition in  a field  of  4700  oersteds,  uirected  in  the  (100)  plane  at  angles  of  0°,  40°. 

45  , 55  , and  90  to  the  [ooQ  direction.  The  warming  curves  (magnetic  moment 
as  a function  of  temperature  at  constant  applied  field)  are  shown  in  Figs.  3 to  5. 
When  the  field  is  within  40°  of  the  [00l]  direction,  this  [oOl]  becomes  the  c axis. 


Magnetization  (arbitrary  units) 
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Field  (oersteds) 

Fig.  2.  Magnetization  curves  in  the  cubic  [ooij  direction  at  78  K 
after  cooling  with  the  indicated  fields  parallel  to  [j)01j  . 
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When  the  field  is  within  40°  of  the  [do]  direction,  i.  e. . more  than  50  from 
Cool]  , £oio]  becomes  the  c axis  and  the  Cool]  direction  is  magnetically  hard. 
When  the  cooling  field  is  applied  at  45°  to  [oOl]  , the  warming  curves  show 
that  in  approximately  one  half  the  sample  [oOl]  is  the  c axis  and  in  the  re- 
mainder [oio]  is  the  c axis. 

The  effect  of  applying  a magnetic  field  at  various  angles  in  the  cubic 
(110)  plane  was  investigated,  using  the  quasi-static  B-H  loop  tracer  (see  below) 
and  an  oblate  spheroidal  sample  prepared  from  a (110)  slice  of  a magnetite  crys- 
tal. Magnetization  curves  were  measured  in  the  [001]  direction  of  this  sample, 
after  cooling  in  a field  of  10,  000  oersteds  directed  in  the  (110)  plane  at  angles 
of  0°,  25°,  and  40°  to  QoOl]  , The  identity  of  the  curves  shows  that,  in  each 
case,  the  c axis  is  established  along  the  Tool]  direction.  All  these  results 
are  consistent  with  the  picture  that  the  c axis  is  established  along  the  cube  edge 
nearest  the  direction  of  the  magnetization  in  each  domain.  Generally  similar 
results  have  been  obtained  recently  with  natural  crystals  of  magnetite  by 
Williams  and  Bozorth.  ^ 


SWITCHING  OF  THE  MAGNETIC  AXIS 
The  phenomenon  of  axis  switching,  i.e.,  the  change  of  the  c axis  from 
one  cube  edge  to  another  below  the  transition,  was  first  observed  by  Bickford  ^ 
in  the  microwave  resonance  experiments  at  liquid  nitrogen  temperature.  The 


effect  is  apparent  in  the  warming  curves  show:,  in  Fig.  5.  The  magnetic  axis 
was  established  along  -he  [oiOj  direction  so  that  L 0 01 J was  magnetically  hard, 


as  shown  by  its  small  magnetic  moment.  Then,  at  some  temperature  below 


the  transition,  there  is  an  abrupt  increase  in  the  magnetic  moment  indicating 
that  the  [00l]  direction  has  suddenly  become  an  easy  direction.  This  effect  is 
13)  H.  J.  Williams  and  R.  M.  Bozorth,  Revs.  Mod.  Phys.  25,  79  (1953). 
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of  considerable  importance  for  two  reasons.  Obviously  it  must  be  considered 
in  interpreting  experimental  results  and  also  in  planning  experiments,  since 
some  measurements  which  appear  to  be  straightforward  are  impossible  except 
at  very  low  temperatures.  It  also  provides  an  additional  means  of  obtaining 
information  about  the  energies  involved  in  the  transition.  The  temperature  at 
which  this  change  occurs  depends  on  the  external  field  applied  parallel  to  the 
Q 0 OlJ  direction.  The  external  field  necessary  to  cause  axis -switching  at 
various  temperatures  was  determined  by  cooling  the  (100)  spheroid  through  the 
transition  with  a field  of  4700  oersteds  along  the  [oio]  direction  and  then 
measuring  the  warming  curves  with  applied  fields  of  from  1000  to  5000  oersteds. 
A number  of  these  curves  are  shown  in  Fig.  6. 

Axis  switching  was  also  obser/ed  in  certain  directions  in  the  (110)  plane. 
Figure  7 shows  the  magnetization  curves  observed  in  the  [no]  direction  after 
cooling  in  a field  of  10,  000  oersteds  parallel  or  perpendicular  (in  the  (HO)  plane) 
to  this  direction.  An  unexpected  result  is  that  axis-switching  occurred  only  in 
the[l!o]dir ection  and  in  a direction  at  15°  to  it.  From  measurements  of  the 
change  in  length  of  a [100]  bar  (see  Appendix  2)  accompanying  the  axis  switch- 
ing, it  was  evident  that  relatively  large  changes  in  dimensions  (b£/£  ^10  ^)  and 
hence  in  elastic  energy  are  involved  in  this  process.  This  elastic  energy 
apparently  suffices,  in  the  case  just  mentioned,  to  prevent  the  axis-switching 
from  occuring  at  the  field  indicated  by  Eq.  (2),  although  it  may  still  occur  in 
a sufficiently  high  field. 

Transfer  of  the  c axis  from  one  cube  edge  to  another  requires  that  elec- 
trons jump  from  ferrous  to  neighboring  ferric  ions  in  the  tetrahedra  formed  by 
the  iron  ions  in  the  octahedral  interstices  (cf.  Fig.  1).  A crystal  in  which  the 
c axis  has  been  established  along  the  [oio]  direction,  placed  in  a magnetic 
field  directed  along  the  [oOl]  direction,  could  reduce  its  energy  considerably 


if  [OOl}  became  the  easy  direction.  The  dependence  of  this  energy  difference 
on  the  magnetic  induction  is  complicated  since  it  is  determined  by  the  aniso- 
tropy energy. 

At  present,  we  will  be  content  with 


a statistical  approach  to  this  problem 
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We  consider  an  electron  and  its  two  possible  positions,  an  "a"  site  correspond- 
ing to  the  c axis  in  the  [oio]  direction  and  a "b"  site  corresponding  to  the 
c axis  in  the  Q)01^  dlrecli°n*  Initially,  there  is  a high  probability  that  an  elec- 
tron will  be  at  an  "a"  site,  since  the  process  of  cooling  the  crystal  in  a field 
directed  along  [oio]  slightly  lowered  the  energy  of  the  "a"  site  relative  to  that 
of  the  "b"  site.  The  application  of  a magnetic  field  in  the  L001]  direction  lowers 
the  energy  of  site  "b".  We  will  assume  that  the  change  in  the  energy  of  the  "b" 
sites  is  proportional  to  the  internal  magnetic  induction.  The  net  probability  that 

an  electron  will  jump  from  an  "a"  to  a "b"  site  is  then  given  by  the  expression 

pB. 

exp  (-u/kT)  (1  - exp  Q-pB./kTj)  = pjA  exp  (-u/kT)  , (1) 

where  pB^  = the  energy  difference  between  sites  "a"  and  "b"  due  to  the  magnetic 
field.  In  writing  this  expression,  we  have  assumed  that  the  initial  energy  dif- 
ference between  the  "a"  and  "b"  sites  is  small  compared  to  pB..  We  have  so 
far  neglected  the  effects  of  the  local  disordering  which  accompanies  a jump 
from  an  ua"  to  a "b"  site.  This  disordering  would  raise  the  energy  of  an  elec- 
tron in  the  "b"  site  and  hence  we  may  expect  an  electron  would  return  to  an  "a" 
site  unless  some  of  its  neighbors  also  jumped  to  "b"  sites  before  it  had  time  to 
return.  We  may  say  crudely  that  the  net  probability  (Eq.l)  must  be  equal  to 
some  fixed  value  before  the  "b"  sites  become  stable  and  the  axis  switches. 

Thus  we  can  write  the  relation  between  the  magnetic  induction  and  the  temper- 
ature at  which  the  c axis  switches  as 

B.  = CT(exp  u/kT)  (2) 

where  B.^  is  the  internal  induction  prior  to  the  axis -switching  and  T is  the  temper- 
ature at  the  mid-point  of  the  swishing.  Figure  8 shows  the  experimental  data 
calculated  from  the  curves  shown  in  Fig.  6.  It  is  clear  that  the  data  fits  the 
above  equation  with  u = 0.  033  cv. 
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Fig.  8.  Induction  necessary  to  switch  the  c axis. 


CONDUCTIVITY  OF  MAGNETITE  AND  ITS  RELATION 
TO  THE  TRANSITION 

If  the  conductivity  of  magnetite  is  due  to  the  exchange  of  electrons  between 
the  ferrous  and  ferric  ions  in  the  octahedral  lattice  sites,  we  would  expect  the 
random  distribution  above  the  transition  to  cause  isotropic  conductivity  but  below 
the  transition  the  conductivity  should  be  anisotropic  (cf.  Fig.  13).  This  aniso- 
tropy cod  its  dependence  on  both  crystallographic  orientation  and  the  direction 
of  a magnetic  field  during  cooling  have  been  studied  in  detail.  The  conductivity 
was  measured  by  the  usual  four -terminal  method.  The  sample  holder  (Fig.  9) 
was  designed  to  maintain  the  sample  and  the  contacts  in  as  nearly  as  possible  an 
isothermal  region  in  order  to  minimize  spurious  thermal  emfs.  The  temperature 
was  measured  using  a copper -constantan  thermocouple  and  a potentiometer  similar 
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Fig.  9.  Conductivity  sample  holder. 
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to  one  described  by  Teele  and  Schuhmann.  ' The  sample  holder  was  mounted 
in  an  evacuated  glass  tube.  Helium  gas,  admitted  through  a needle  valve, 
allowed  the  pressure  to  be  varied  from  10  ^ to  10  ^ mm  and  consequently  the 
warming  and  cooling  rates  could  be  controlled  from  3°  to  50°C  per  hour.  The 
samples  used  were  cut  from  two  synthetic  crystals  in  the  form  of  rectangular 
bars  about  1 cm  in  length  and  with  a cross-sectional  area  of  about  3 sq.  mm. 

The  length  of  the  bar  was  parallel  to  the  cubic  [100],  [no]  , and  [m]  di  - 
rections,  respectively. 

The  conductivity  of  two  of  the  specimens  between  room  temperature  and 

the  transition  is  shown  in  Fig.  10.  The  data,  taken  during  warming  as  well  as 

cooling  cycles,  coincide;  there  was  no  thermal  hysteresis  connected  with  the 

transition.  The  conductivity  is  250  ohm  * cm  ^ at  room  temperature  and  has  a 

very  broad  maximum  at  ca.  15°C,  which  does  not  show  clearly  in  the  figure. 

2 4) 

Previous  investigators  ’ ' have  reported  this  maximum  at  a somewhat  higher 
temperature  (ca.  80°C).  The  transition  is  marked  by  a decrease  of  a factor  of 
^0  in  the  conductivity  in  a temperature  interval  of  ca.  1°.  The  midpoint  of  this 
interval  occurred  at  a temperature  of  119. 4 +_  0.  3°K  (-153.  8°C). 

The  conductivity  below  the  transition  for  samples  cooled  in  zero  magnetic 
field  is  isotropic  (Figs.  11  and  12)  but  does  not  follow  a simple  exponential  law 

<r  = A exp  (-u/kT) 

If  fitted  to  this  type  of  law  over  certain  ranges  of  temperature,  one  obtains  the 
values  of  the  constants  A and  u as  shown  in  Table  1. 

When  the  samples  are  cooled  through  the  transition  in  a magnetic  field, 
the  conductivity  is  anisotropic  and  depends  on  the  orientation  of  both  the  sample 
and  the  magnetic  field  (Fig.  13).  The  fields  of  ca.  9000  oersteds  used  were  more 

14)  R.  P.  Teele  and  S.  Schuhmann,  J.  Research  Nat.  Bur.  Standards  22, 

431  (1439). 
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Fig.  11.  Conductivity  below  the  transition  after  cooling  in  zero  rr  agnetic  field. 
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Table  1.  Conductivity  of  magnetite  below  the  transition. 


Temperature 

A(ohm  * cm  *) 

u(ev) 

78  to  90°K 

8200 

0.11 

56  to  77° 

1400 

0.09 

40  to  52° 

1.  5 

0.06 

than  sufficient  to  saturate  the  samples  above  the  transition.  The  anisotropy 
increases  rapidly  from  the  transition  to  approximately  95°K,  and  then  re- 
mains practically  constant  at  lower  temperatuz  es.  This  anisotropy  is  quite 
large,  the  ratio  of  maximum  to  minimum  conductivity  being  as  large  as  1.  88. 
Table  2 gives  the  ratios  of  the  conductivities  for  a number  of  combinations  of 
sample  orientation  and  magnetic  field  direction. 


Table  2.  Ratios  of  conductivities  of  magnetite  in  different 
directions  below  the  transition 


T 

X 

n? 

o 

L^j 

b 

(r[l00]  ||  H 

jE 

o 

o 

b 

<r  [lio]  45°H 

<t  [lOo]  ||  H 

<r[l00]lH 

<r[U0]±H 

<r  [l  lo]  45°H 

<r[ilo]  -LH 

-[Hi] 

1 1 17.  6°K 

1.44 

1.40 

1.25 

1.1  2 

- 

111.1 

1.63 

1.51 

1.  21 

1.  25 

- 

105.  3 

1.79 

1.  59 

1.  27 

1.  25 

- 

100.0 

1.87 

1.39 

1.  24 

1. 28 

- 

95.  2 

1.88 

1.61 

1.25 

1.  28 

1.  56 

90.9 

- 

1.61 

1.  27 

1.  27 

- 

87.0 

- 

1.60 

1.  26 

1.  27 

- 

83.  3 

- 

1.60 

1.  26 

1.  26 

The  observed  anisotropy  in  the  conductivity  agrees  well  with  that  which 
can  be  predicted  from  Verwey's  model  of  the  lew  temperature  structure  of 


magnetite.  Although  the  exchange  of  electrons  between  ferrous  and  ferric  ions 
has  become  much  more  difficult  below  the  transition,  it  is  still  the  dominant 
mechanism  of  conduction.  In  the  ordered  structure  this  will  lead  to  an  aniso- 
tropy since  ferrous  and  ferric  ions  occupy  adjacent  sites  in  only  four  of  the 
cubic  [lio]  dir  ections.  Although  the  ordered  structure  has  orthorhombic  sym- 
metry, its  conductivity  anisotropy  corresponds  to  tetragonal  symmetry  (neg- 
lecting the  small  departure  of  the  b/a  axial  ratio  from  unity).  The  conductivi- 
ty of  a tetragonal  crystal  in  an  arbitrary  direction  is  given  by 

cr  = o-jj  sin^O  + o-jj  cos^  0 , (3) 

where  9 is  the  angle  between  the  direction  and  the  c axis,.  ^e  conductivity 
along  the  c axis,  and  <r^  the  conductivity  normal  to  the  c axis.  It  can  be  seen 
easily  from  the  ordered  structure  that  the  ratio  = 2.  Thus  we  can  write 

the  conductivity  below  the  transition  as 

«r  = B (1  + cos^  0)  (4) 

In  case  the  order  is  not  perfect,  we  would  expect  the  conductivity  to  contain,  in 
addition,  an  isotropic  term 

o-  = A + B (1  + cos^  0)  . (5) 

The  quantity  B/(A  + B)  should  be  a measure  of  the  amount  of  order.  The  values 
of  A,  B,  and  B/(A  + B)  derived  from  our  conductivity  data  for  a number  of  temper 
atures  are  given  in  Table  3. 


Table  3.  Anisotropy  of  conductivity  below  the  transition. 


T 

A 

B 

B(A  + B) 

1 17.  6°K 

0. 2660 

0.1  270 

0.324 

111.1 

0.0718 

0. 0741 

0.  508 

10*.  3 

0. 0202 

0.0396 

0.  662 

1 00.  0 

0. 0079 

0.0194 

0.711 

95.  2 

0. 0026 

0.  01155 

0.81  6 
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This  table  shows  that  with  decreasing  temperature  the  order  increases, 
as  one  would  expect,  since  the  ordering  process  involves  an  activation  energy. 


MAGNETIC  ANISOTROPY 

The  magnetic  anisotropy  energy  is  the  difference  between  the  energy 
required  to  magnetize  a sample  in  an  arbitrary  direction  and  in  the  direction 
of  easy  magnetization.  The  angular  dependence  of  this  anisotropy  energy  is 
determined  by  the  symmetry  of  the  crystal.  In  the  ordered  structure  of  mag- 
netite a determination  of  the  angular  dependence  of  the  anisotropy  energy  is 
particularly  important  since  the  actual  deformation  of  the  crystal  structure  at 
the  transition  is  very  small  and  therefore  difficult  to  determine  directly,  while 
the  anisotropy  energy  below  the  transition  is  very  large.  Consequently,  the  de- 
termination of  the  anisotropy  energy  provides  a sensitive  method  for  checking 
the  crystal  symmetry. 

The  method  of  calculating  the  dependence  of  the  anisotropy  energy  on 
orientation  for  any  type  of  crystal  symmetry  is  well  known.  The  anisotropy 
energy  is  expanded  as  a power  series  in  the  direction  cosines  relative  to  the 
crystallographic  axes,  only  those  terms  consistent  with  the  symmetry  being 
retained.  The  existence  of  a plane  of  symmetry  normal  to  a crystallographic 
axis  requires  the  direction  cosine  relative  to  that  axis  to  appear  only  in  even 
powers.  Since  the  ordered  structure  of  magnetite  belongs  to  the  orthorhombic 
holohedral  crystal  class,  it  has  three  such  symmetry  planes  (one  normal  to 
each  axis).  The  remaining  symmetry  elements  of  this  crystal  class,  i.e.,  the 
three  two-fold  axes  and  the  center  of  symmetry,  can  be  compounded  from 


appropriate  sequences  of  reflections  in  the  three  symmetry  planes  and  there- 
fore do  not  impose  any  additional  restrictions  on  the  form  of  the  anisotropy 
energy.  Thus  we  have 


(6) 
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where  the  K's  are  the  anisotropy  constants,  and  the  a's  the  direction  cosines 
relative  to  the  crystallographic  axes.  The  identity  + a*  + = 1 has  been 

used  to  eliminate  all  the  terms  containing  a^.  For  magnetite  below  the  trans- 
ition, me  c axis  is  the  easy  direction  and  one  of  the  a and  b axes  is  much 
harder  than  the  other.  Since  this  difficult  axis  could  not  be  experimentally 
identified  with  one  of  the  crystallographic  axes,  the  b axis  was  arbitrary 
chosen  as  the  difficult  one.  This  choice  may  be  inconsistent  with  the  assign- 
ment of  the  true  crystallographic  axes  but,  if  this  should  prove  to  be  the  case, 
only  the  subscripts  1 and  2 of  the  anisotropy  constants  would  have  to  be  inter- 
changed. 

The  anisotropy  constants  were  computed  from  the  areas  between  the 

magnetization  curves  in  different  crystallographic  directions.  Because  of  the 

high  fields  needed  to  saturate  magnetite  below  the  transition,  it  was  necessary 

to  obtain  the  magnetization  curves  in  fields  as  high  as  15,  000  oersteds.  Since 

the  pendulum  magnetometer  is,  at  present,  limited  to  fields  of  less  than  5000 

15) 

oersteds,  the  quasistatic  B-H  loop  tracer  constructed  in  this  laboratory  ' 
was  modified  to  trace  the  magnetization  curves  of  samples  mounted  in  the  gap 
cf  the  large  electromagnet.  The  sample  was  mounted  at  the  center  of  a fixed 
coil.  The  contribution  of  the  field  itself  to  the  flux  linking  this  coil  was  elimi- 
nated by  mounting  a second  coil  in  the  gap,  some  distance  from  the  sample 
(Fig.  14).  The  two  coils  were  connected  in  opposition  and  were  adjusted  to 
give  as  nearly  as  possible  zero  output  with  no  sample  in  position.  The  voltage 
of  these  two  coils  was  integrated  and  applied  to  one  axis  of  a recorder.  The 
second  axis  was  driven  by  the  output  of  a rotating  coil  gaussmeter,  also 
mounted  in  the  magnet  gap.  The  curves  obtained  on  the  recorder  are  only 
proportional  (not  equal)  to  the  magnetization,  because  of  the  flux  leakage  through 
15)  D.  J.  Epstein  and  3.  Frackiewicz,  Progress  Rep.  XI,  May,  1952,  p.  55. 
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Fig.  14.  Empirical  relation  for  the  conductivity  above  transition. 
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the  small  air  gap  between  the  sample  and  the  coil.  Since  the  saturation  mag- 
netization of  magnetite  at  room  temperature  has  been  measured  by  a number  • 
of  investigators,  this  was  used  to  calibrate  the  scale.  The  sample  used  for 
these  measurements  was  an  oblate  spheroid  cut  from  a (110)  plane  slice  of  a 
synthetic  crystal  of  magnetite.  It  was  shaped  by  a method  developed  in  this 
laboratory  so  that  it  was  a much  better  approximation  to  an  oblate  spheroid 
than  the  (100)  sample.  Its  axes  were  3.84  and  0.97mm,  respectively,  and  the 
demagnetizing  factor  calculated  from  these  dimensions  was  1.  82  + . 08,  as  com- 
pared to  1.75  computed  from  the  [in]  magnetization  curve  at  room  temper - 
acure. 

Room  temperature  magnetization  curves  were  measured  for  several 

2 5 

reasons:  the  saturation  magnetization  at  room  temperature  (477  cgs  units/cc  ’ ) 
was  used  to  fix  the  scale  of  the  magnetization;  these  curves  provided  a check 
on  the  demagnetizing  factor  calculated  from  the  dimension  of  the  sample;  and, 
finally,  a comparison  of  the  anisotropy  constant  at  room  temperature  with  the 
values  obtained  by  other  methods  provided  an  experimental  check  on  the  accu- 
racy and  reliability  of  the  equipment.  For  cubic  holohedral  crystals,  the 

anisotropy  energy  can  be  written  as 

_ v , 2 2 22  2 2 » v 222 

Ea  = Kl(ala2  + a2a3  + a3al)  + K2alQ2a3  * (7) 

The  value  of  at  room  temperature  was  -1.12  + 0.  05  x 10^  ergs/cc,  as  de- 
termined from  the  areas  between  the  QOO]  and  [lio]  as  well  as  between  the 
f 100 J and  Tllll  curves.  This  is  in  good  agreement  with  the  value  of  -1.12 

for  from  microwave  resonance  measurements  ^ and  -1.  22  from  torque  mag- 

13) 

netometer  measurements.  ' is  zero,  within  the  limits  of  experimental  error. 

The  low  temperature  curves  were  obtained  by  immersing  the  sample  and 
the  surrounding  coil  in  liquid  air.  The  sample  reached  equilibrium  very  quick- 
ly at  a temperature  of  ca.  -188°C.  The  saturation  magnetization  of  magnetite 
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at  this  temperature  was  508  cgs  units/ cc,  in  excellent  agreement  with  the 
value  of  510  at  -195°  obtained  by  Domenicali.  ^ Magnetization  curves  were 
obtained  in  the  cubic  [lOO]  . [lio]  , and  [ill]  directions  and  also  15°,  30°, 

45°  and  75°  to  [lOoj  after  cooling  with  either  a parallel  or  a perpendicular 
magnetic  field  of  ca.  10,  000  oersteds.  Curves  were  also  obtained  in  the 
cubic  [lOo]  and  [lioj  directions  after  cooling  in  zero  field. 

The  twinning  that  occurs  in  magnetite  when  it  is  cooled  through  the 
transition  introduces  a number  of  complications.  A magnetic  field  applied  to 
the  sample  as  it  is  cooled  will  align  the  c axis  along  the  cube  edge  nearest  the 
applied  field.  Since  the  anisotropy  energy  is  considerably  lower  along  the 
a axis  than  along  b,  the  a axis  will  tend  to  lie  closer  to  the  external  field.  ^ 
These  two  conditions  are  sufficient  to  remove  the  twinning  only  in  certain 
orientations.  In  some  cases,  particularly  when  the  applied  field  makes  nearly 
equal  angles  with  two  or  three  cube  edges,  there  is  uncertainty  about  the  re- 
lative volumes  of  the  sample  which  will  assume  each  of  the  possible  orienta- 
tions. A much  more  serious  difficulty  connected  with  this  twinning  is  the 
interactions  that  occur  between  the  different  regions  of  the  sample.  Thus,  in 
many  cases,  a twinned  sample  cannot  be  treated  by  a simple  superposition  of 
the  various  orientations  which  it  contains.  These  interactions  are  discussed 
below  in  detail. 

The  anisotropy  constants  were  calculated  by  trial  and  error  from  the 
areas  between  the  magnetization  curves.  Due  to  the  uncertainties  introduced 
by  the  possible  twinning  and  the  interactions  between  the  different  regions  of 
the  sample,  a simple  fitting  of  the  measured  areas  between  the  various  mag- 
netization curves  did  not  appear  to  be  a sufficient  check  on  the  values  of  the 

16}  H.  J.  Williams,  R.  M.  Bozorth  and  M.  Goertz,  "Mechanism  of  the  Trans- 
ition in  Magnetite  at  Low  Temperatures,  " Phys.  Rev.  , in  press. 
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anisotropy  constants.  The  constants  obtained  from  these  areas  were  there- 
fore used  to  compute  the  magnetization  curves  for  the  cubic  [lOo]  and  fill! 
directions  after  cooling  in  a perpendicular  field.  The  constants  were  then 
adjusted  to  improve  the  fit  between  the  computed  and  the  experimental  curves. 
The  values  of  the  constants  finally  obtained  were:  =4.0,  = 9.0,  = 2.2, 

Kj2  = 13.0,  K22  = 10.2  x 105  ergs/cc.  The  observed  and  calculated  areas  are 
compared  in  Table  4,  which  also  shows  the  orientations  of  the  orthorhombic 
axes  used  in  the  computations.  The  calculated  and  experimental  curves  for 
the  cubic  £l00]  ^d  |_illj  directions  are  shown  in  Figs.  15  and  16.  The  aniso- 
tropy energy  is  plotted  as  a function  of  crystallographic  orientation  in  Fig.  17. 

Williams  and  Bozorth^’  have  measured  torque  curves  on  samples  of 
natural  magnetite  at  liquid  nitrogen  temperature,  i.  e. , at  about  -195°C.  When 
their  results  are  expressed  in  terms  of  our  constants,  as  shown  in  Table  5,  the 
agreement  between  the  two  sets  of  constants  is  relatively  good,  considering  the 
differences  in  the  temperature  and  the  crystals  used.  Generally,  the  natural 
crystals  are  not  exactly  stoichiometric  and  consequently  would  not  have  as  per- 
fect an  ordered  structure  as  the  synthetic  crystals.  The  amount  of  order  in- 
creases rapidly  as  a crystal  is  cooled  below  the  transition.  Thus  the  7° differ- 
ence in  temperature  between  the  two  sets  of  measurements  partially  compensates 


Table  5.  Anisotropy  constants  of  magnetite  at  low  temperature. 


Constants 

/ -5 

ergs/cc  x 10 

Method 

Magnetization  curves 

Torque  curves 

*1 

4.  0 

2.  8 

K2 

9.0 

7.5 

*11 

2.  2 

2.  2 

*12 

13.  0 

12.4 

*22 

10.  2 

7.4 
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Fig.  17.  Anisotropy  energy  of  magnetite,  showing  the  path  followed 
8 by  the  magnetisation  when  a field  is  applied  in  the  ortho 

rhombic  \U0\  direction.  The  projections  of  the  ortho- 
rhombic  a,  b,  and  c axes  are  at  the  corners  (a),  (b),  and  (c) 
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for  the  inherent  difference  in  the  crystals.  The  deviations  in  all  constants 
are  in  the  same  direction.  Their  effect  on  the  magnetization  and  torque 
curves  is  cumulative;  obviou,ly  there  is  a real  difference  between  the  natural 

and  the  synthetic  crystals. 

INTERACTIONS  BETWEEN  THE  REGIONS 
OF  A TWINNED  CRYSTAL 

When  a sample  of  magnetite  is  cooled  through  the  transition  in  such  a 
way  that  twinning  occurs,  i.e.,  the  orthorhombic  axes  do  not  have  the  same 
orientation  throughout  the  sample,  we  have  a very  complicated  problem  to 
deal  with.  The  simplest  assumption  that  can  be  made  is  that  the  regions  of 
the  sample  in  which  the  axes  are  oriented  differently  act  completely  indepen- 
dently of  each  other.  Thus  the  behavior  of  the  whole  sample  would  be  a simple 
superposition  of  the  behavior  of  the  individual  regions.  Such  a simplified  pic- 
ture ignores  the  interactions  between  the  different  regions,  which  could  arise 
from  a number  of  sources;  the  magnetizations  in  neighboring  regions  exert 
torques  on  each  other  which  would  affect  the  motions  under  an  external  field; 
each  region  would  be  strained  owing  to  the  distortions  experienced  by  its 
neighbors;  and  free  poles  on  the  boundaries  between  the  regions  could  cause 
discontinuities  in  the  internal  field.  The  effects  due  to  the  last  two  causes 
could  be  very  complex,  and  lacking  the  information  needed  to  even  estimate 
them,  we  can  only  assume  that  they  do  not  greatly  alter  the  magnetization 
curves. 

The  effect  of  torques  exerted  by  the  magnetization  in  one  region  on  that 
in  neighboring  regions  is  particularly  large  in  two  of  the  experimental  mag- 
netization curves:  the  cubic  [lOO]  direction  cooled  in  zero  field,  and  the  cu- 
bic [ill]  direction  cooled  in  a parallel  field.  When  a sample  is  cooled  through 
the  transition  in  demagnetized  state,  the  cubic  [lOO]  direction  will  become  the 
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orthorhombic  [ool]  in  l/3  of  the  sample;  the  orthorhombic  [lio]  in  l/3;  and 
the  orthorhombic  [lio]  in  l/3,  since  the  orthorhombic  a and  b axes  develop 
along  face  diagonals  of  the  cubic  unit  cell.  The  magnetization  curves  in  the 
orthorhombic  [lio]  and  Qllo]  directions  will  be  identical.  The  orthorhombic 
[lio]  curve  can  be  obtained  by  cooling  with  a field  normal  to  the  cubic  [lOO] 
direction,  and  the  [oOl]  curve  with  a parallel  field.  Since  the  c axis  (001)  is 
the  easy  direction,  if  there  were  no  interaction,  the  magnetization  curve  for 
the  cubic  [ lOo]  direction,  after  cooling  in  zero  field,  should  be  given  by 

h0  = i io  = yj  * aji  • <*> 

All  these  quantities  can  be  obtained  from  the  curves  in  Fig.  15.  The  curve 
calculated  from  this  equation  is  compared  in  Fig.  18  with  the  experimental 
curve.  The  discrepancy  is  due  to  the  interaction  between  the  different  regions, 
which  causes  an  additional  torque  to  affect  the  magnetization.  In  part  of  the 
sample,  there  is  magnetization  parallel  to  the  field  and  in  the  remainder  of 
the  sample,  a magnetization  I,  at  an  angle  0 to  the  field.  The  magnetization 
parallel  to  the  field  will  exert  a torque  on  the  remainder  of  the  magnetization. 
The  energy  arising  from  this  torque  will  be  of  the  form  l\  cos  0.  Thus  the 
energy  of  a region  whose  c axis  is  perpendicular  to  the  applied  field  will  be 

E = E - HI  cos  0 - Cl*  cos  0 . (9) 

X as  3 

where  Ea  i.  the  anisotropy  energy.  H the  internal  field,  1,  the  .atur.tion  mag- 
netisation, and  C •«  empirical  con.tant.  To  find  the  equilibrium  position  of  the 
magnetization,  we  set  dE d0  = 0 and  obtain 

H - Cl  = -dE  / d0/l  »in  0 • (10) 

s a s 

The  corresponding  equation  for  a sample  containing  only  regions  with  their 
c axes  normal  to  the  applied  field  is  the  same  except  that  the  term  -CIs  is  not 
present.  Therefore,  the  magnetization  curve  for  the  cubic  [loo]  direction 


- 36  - 


- 37  - 


cooled  in  zero  field  will  be  given  by  the  equation 

H0  = hj.  - ci,  : io=  V3  + 2Ix/3  (n) 

The  fit  between  this  equation  with  C * 1.  and  the  experimental  curve  is  quite 
satisfactory  (cf.  Fig.  18).  The  discrepancy  in  fields  below  800  oersteds  is 
due  to  strains  induced  in  the  sample  as  it  cooled  through  the  transition,  and 
to  the  effects  of  the  magnetostriction.  A comparison  of  this  curve  with  that 
for  the  cubic  ^100]  direction  cooled  in  a parallel  field  indicates  that  these 
effects  are  much  larger  here,  as  we  would  expect. 

This  case  is  one  in  which  the  interaction  has  a very  large  effect  on  the 
magnetization  curve  for  two  reasons:  the  interaction  torque  is  large  because 
the  angle  between  the  magnetizations  is  large;  and  this  torque  acts  in  the  same 
direction  as  the  torque  due  to  the  applied  field.  Another  case  in  which  these 
conditions  are  present  is  that  of  the  cubic  direction  cooled  in  a parallel 

field.  In  this  case,  the  c axis  will  be  distributed  over  all  three  cube  edges 
with  tbe  a and  b axes  so  oriented  that  the  cooling  field  lies  in  an  a-c  plane 
throughout  the  sample.  The  interaction  is  more  complicated  here  because  each 
magnetization  will  be  subject  to  two  interaction  torques  and,  in  addition,  the 
angle  between  two  magnetizations  is  not  the  same  as  the  angle  between  each 
one  and  the  applied  field.  We  can  see  the  effect  of  the  interaction  directly  by 
comparing  the  cubic  [mj  curves  after  cooling  in  different  fields.  After  cool- 
ing in  a perpendicular  field,  the  cubic  (110)  plane  will  transform  into  the  a-c 
plane  of  the  orthorhombic  crystal,  with  no  multiple  orientation.  Thus  the  dif- 
ference in  the  two  cubic  [m]  curves  (Fig.  16)  is  due  to  the  interaction  occur 
ring  after  cooling  in  a parallel  field.  In  most  of  the  other  cases  we  have  ob- 
served, the  effects  of  the  interactions  will  be  much  smaller. 

These  interactions  probably  account  for  the  failure  of  the  anisotropy 
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Fig.  18.  Magnetization  curve  at  85°K  in  the  £l00}  ^rec^on  after  cooling  in  zero  field. 
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constant*  obtained  by  static  methods,  to  explain  the  microwave  resonances 
observed  by  Bickford.1)  The  twinning  that  occurs  at  the  transition  accounts 
qualitatively  for  the  multiple  resonances  observed,  except  in  the  case  of  the 
cubic  (100)  plane.  Some  of  the  resonances  in  this  plane  may  be  due  to  lack  of 
magnetic  saturation  in  the  sample.  When  the  observed  anisotropy  constants 
are  used  to  calculate  the  positions  of  the  resonance  peaks,  these  calculated 
values  do  not  agree  with  those  observed.  Interactions  of  the  type  described 
above  produce  torques  which  would  affect  the  precession  of  the  magnetization 
and  hence  the  position  of  the  resonance.  It  will  be  necessary  to  perform 
these  experiments  again  under  conditions  such  that  no  twinning  occurs  in  order 
to  determine  whether  the  static  anisotropy  constants  will  account  for  the  ex- 
perimental results. 

SUMMARY  AND  CONCLUSIONS 

The  orthorhombic  structure  of  magnetite  below  the  transition  is  con- 
firmed by  a wide  variety  of  experimental  evidence.  Direct  determination  of 
the  unit  cell  size,  both  by  X rays  ) and  strain-gauge  measurements, 11)  show 
that  the  orthorhombic  axes  are  oriented  so  that  the  c axis  lies  along  one  of  the 
original  cube  edges  and  the  a and  b axes  along  the  perpendicular  cubic  face- 
diagonal  directions.  The  magnetic  anisotropy  energy  gives  perhaps  even  more 
convincing  evidence  since  it  is  so  large  that  there  can  be  no  doubt  about  its 
form,  whereas  the  changes  in  the  unit  cell  itself  are  very  small.  The  aniso- 
tropy of  the  conductivity  agrees  with  that  to  be  expected  from  Verwey's  ordered 
structure. 

In  addition  to  accounting  for  the  anisotropic  properties  below  the  trans- 
ition, Verwey's  model  explains  the  influence  of  a magnetic  field,  applied  during 
cooling  on  the  properties  of  the  crystal  below  the  transition.  The  axis -switch- 
ing effect  was  not  predicted  from  this  model,  but  the  model  does  provide  an 


explanation  for  it.  While  it  is  true  that  the  ordering  has  not  been  directly 
observed,  e.g. , superlattice  lines  have  not  yet  been  found  in  X-ray  or  neu- 
tron diffraction  experiments,  it  accounts  so  completely  and  quantitatively  for 
the  behavior  of  magnetite  at  and  b low  its  low-temperature  transition  that 
there  can  be  little  doubt  of  the  correctness  of  Verwey's  hypothesis. 

The  twinning  which  occurs  at  the  transition  greatly  complicates  the 
magnetic  behavior  at  lower  temperatures  because  the  different  regions  of  a 
twinned  crystal  interact  in  a complex  way  with  each  other.  These  inter- 
actions are  quite  large  in  some  cases  and  probably  account  for  the  failure  of 
the  static  anisotropy  constants  to  explain  the  microwave  resonance  data. 

One  remaining  problem  of  major  interest  concerning  the  low-temper- 
ature region  of  magnetite  is  a detailed  interpretation  of  the  conductivity  data. 
The  increase  in  conductivity  anisotropy  with  decreasing  temperature  (Table  3) 
appears  to  be  directly  related  to  the  increase  in  long-range  order.  A greatly 
over-simplified  band  model,  which  exhibits  a discontinuity  in  the  conducti- 
vity such  as  that  in  magnetite  at  the  transition,  is  described  in  Appendix  3. 
Such  models  must  necessarily  be  very  speculative  until  much  more  experi- 
mental data  become  available  particularly  on  the  Hall  effect  and  the  thermo- 
electromotive force. 
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Appendix  1 . 

DETAILS  OF  THE  CALCULATION  OF  THE  LOW  TEMPERATURE 

MAGNETIZATION  CURVES 

The  calculation  of  the  magnetization  curve  in  the  cubic  [m]  direction 
of  a sample  cooled  in  a magnetic  field  normal  to  this  direction  in  a (110)  plane, 
proceeds  as  for  a cubic  crystal.  Both  the  easy  (c-axis)  direction  and  the 
applied  field  lie  in  the  a-c  plane  of  the  orthorhombic  structure  and  there  is  no 
twinning  to  complicate  the  problem.  Since  the  anisotropy  energy  is  symmetric 
about  this  plane,  there  will  be  no  torque  acting  to  move  the  magnetization  out 
of  the  plane  and  its  equilibrium  position  will  be  determined  by  the  condition 
that  the  total  energy 

ET  = Ea  - HIs  cos  (54°44'  - 8)  (12) 

is  a minimum  (9  i.i  the  angle  between  1 and  c axis).  The  anisotropy  is  given 
by 

Ea  = sin20  + *11  sin40  . (13) 

Setting  dE.j./d0  = 0,  we  obtain 

H = (Kj  sin  20  + 2K][1  sin20  sin  20)/lg  sin(54°44'  - 0)  , (14) 

and  we  have 

1 = cos  (54°44'  - 0)  . (15) 

From  these  parametric  equations,  the  magnetization  curve  can  be  computed. 

The  calculation  for  a cubic  [lOo]  direction  cooled  with  a field  along  a 
perpendicular  [lio]  is  much  more  complicated.  We  will  neglect  the  interactions 
between  the  different  regions,  since  the  experimented  curves  (Fig.  15) 


- 41  - 


indicate  that  in  this  case  they  do  not 
significantly  affect  the  magnetization 
curve.  This  curve  will  be  that  for  the 
orthorhombic  [lio]  direction.  The 
anisotropy  energy  is  not  symmetric 
about  the  orthorhombic  (110)  plane, 
and  consequently  there  will  be  a torque 
acting  to  turn  the  magnetization  out  of 
this  plane.  To  find  the  path  followed  by 
the  magnetization,  we  use  the  spherical 
co-ordinates  indicated  in  Fig.  19  and 
express  the  anisotropy  energy  in  terms 


Fig.  19.  Co-ordinate  system  for 
calculation  of  ortho- 
rhombic 110  magneti- 
zation curve. 


af  the  angles  0 and  4>.  Since  the  field 


exerts  no  torque  in  the  direction  of  4>.  the  path  will  be  determined  by  the 


condition 


In  this  case,  we  have 


9 E /3<j>  = 0 

3i 


= ^^-(cos  0 + sin  4>  sin  0) 

(2)  ' (17) 

o = Ly(cos  0 • sin  <J>  sin  0)  , 

2 (2)1/2 

and,  substituting  these  expressions  into  the  anisotropy  energy  (Eq.  6)  and  dif- 


ferentiating, we  obtain 

cos  $ ^sin  [(^1  + *12  + K22^  Sln 
+ sin2  4*  ^3  (K^  - K22)  sin30  cos  0] 

+ sin  4»[(3K11  - Kj2  + 3K22)  sin2  0 cos20  + (Kj  + K2)  sin20j 
♦[(Ku-km>  sin  8 cos3  0 + (Kj  - K2)  sin  0 cos  e3}  = 0.  (18) 

If  4>  / 90°,  this  is  a cubic  equation  in  sin  <J>.  For  values  of  the  anisotropy  con- 
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st&nts  in  the  range  of  those  of  magnetite,  this  equation  has  only  one  real  root. 
We  can  find  any  number  of  points  along  the  path  of  the  magnetization  by  sub- 
stituting various  values  of  0 into  this  equation  and  solving  numerically  for  the 
corresponding  values  of  0.  The  path  is  indicated  in  Fig.  17.  This  method  is 
applicable  for  any  orientation  of  the  applied  field  relative  to  the  crystallo- 
graphic axes  and  for  any  type  of  crystal  symmetry.  In  many  cases,  however, 
a graphical  method  described  by  Bozorth  and  Williams  ^ would  be  less  la- 
borious. 


With  the  path  determined,  we  can  set  up  the  equations  for  the  magneti- 
zation curve.  The  energy  is  given  by 

ET  = Ea  - HIs  cos  0 , (19) 

and,  after  differ entiatin,  we  obtain 


We  then  have 


H = - 3E  /d0/l  sin  0 

fit  8 


I = Ia  cos  0 


The  derivative  3Eft/30  is  given  by 

^Ea/^e  = _1/2  (Kj  + K2)  sin  20  cos2$  + (Kj  - K2)  sin  cos  26 
-i/2  (Kj,  t + K^)  cos20  sin  20 
+l/2  (K^  + sin4<j>  sin20  sin  20 

tl/4  (JKjj  - Kj,  + 3K22>  sin2$  sin  40 

”3/4  (*11  ' Kzz>  8in  * cos24>  sin2  20  + (KJ1  - K22)  sin  * cos40 
- - K22)  sin  4.  sin  0 , (22) 

and  must  be  evaluated  at  points  along  the  path  of  the  magnetization. 
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AposnrJix  2. 

MAGNETOSTRICTION  OF  MAGNETITE  AT  ROOM  TEMPERATURE 

Measurements  of  the  magnetostriction  of  natural  crystals  of  magnetite,1®^ 
indicated  that  its  behavior  was  very  unusual  in  that  the  longitudinal  saturation 
magnetostriction  in  the  [no]  direction  was  not  intermediate  between  that  for 
the  [100]  and  fill}  directions.  The  magnetostriction  of  the  synthetic  crys- 
tals of  magnetite  was  measured  using  the  dilatometer  constructed  by  Domeni- 
cali.  2^  Changes  in  the  length  of  the  sample  alter  the  separation,  and  hence 

the  capacitance,  of  two  metal  plates.  These  changes  in  capacitance  were 

19) 

measured  with  a beat-frequency  equipment  developed  in  this  laboratory.  ' 

With  this  apparatus,  the  minimum  change  in  length  that  could  be  measured 
-7 

was  6 x 10  mm.  The  dilatometer  was  calibrated  by  measuring  the  thermal 
expansion  of  magnetite,  nickel  and  fused  quartz  between  15>J  and  30°C. 

The  samples  used  for  the  magnetostriction  measurements  were  the  rect- 
angular bars  used  for  the  conductivity  measurements.  The  observed  magneto- 
striction depends  in  a complicated  way  on  the  shape  of  the  sample  and  the 
elastic  constants.  2^  In  iron,  this  correction  for  the  "form  effect"  is  extreme- 
ly important  but,  since  it  is  proportional  to  the  square  of  the  magnetization 

-7 

it  is  very  small  in  magnetite  (of  the  order  of  10  ) and  was  ignored.  The  most 

serious  drawback  of  the  bar-.-.haped  samples  is  that  the  demagnetizing  field  is 
not  uniform.  Hence,  the  magnetizafion  processes  are  smeared  out  and  the 
magnetostriction  data  do  not  provide  information  on  the  various  translational 
and  rotational  processes  occurring  during  the  magnetization  of  the  sample. 

Only  the  saturation  magnetostrictions  have  full  significance. 

18)  C.  W.  Heaps,  Phys.  Rev.  24.  60  (1924). 

19)  J.  G.  Jelatis,  Tech.  Rep.  7;  J.  Appl.  Phys.  19,  419  (1948). 

20)  R.  Becker  and  W.  Dttring,  "Ferromagnetismus,  " Springer,  Berlin,  1939,  p.305. 
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The  longitudinal  magnetostriction  was  measured  in  the  three  crystallo- 
graphic directions  [lOo]  , Qio]  , [ill].  The  transverse  magnetostriction  was 
measured  by  rotating  the  dilatometer  90°  in  the  magnet  gap,  so  that  the  field 
was  applied  parallel  to  a short  dimension  of  the  sample  and  the  elongation 
measured  along  the  length.  It  was  apparent  from  the  values  of  the  longitudinal 
saturation  magnetostrictions  that  the  usual  two-constant  expression  for  the  mag- 
netostriction could  cot  be  made  to  fit  the  data.  The  experimental  values  were 
fitted  to  a five-constant  expression  ^ 

bl/£  = - l/3) 

+ h^(2  Oj  a2  P1  p2  + 2 a,  a3  P2  P3  + 2 a3  Uj  P3  P1 ) 

+ h3(s  - 1/3)  + h4(a^p^  + a2p2  + a*  p*  + 2s/3  - l/3) 

+ h5(2a1a2a3pip2  + 2a2a3a^p2p3  + 2a3a1a2p3p1  . (23) 

where  a's  are  the  direction  cosines  ol  the  magnetization,  p's  the  direction  cosines 
of  the  elongation,  and  s=  a ^ t»2  + a2a3  * °3  Q 1 * Since  h3  represents  the  volume 
magnetostriction,  which  is  very  small  '-ompared  to  the  linear  effect,  we  may 
place  h3  = 0. 

The  experimental  results  were  fitted  by  taking 

*>1  = 547>  h2  r m*  h4  = _576»  h5  = 56  x 10'6 
The  calculated  and  observed  values  are  compared  in  Table  6.  The  information 
is  too  scanty  to  permit  these  constants  to  be  very  precisely  determined;  h5, 
particularly,  may  be  in  error  by  a factor  of  2 or  more,  since  it  appears  to  a 
significant  extent  only  in  the  expression  for  the  longitudinal  magnetostriction  in 
the  [ill]  bar. 


21)  R.  Becker  and  W.  Ddring,  op.  cit. , 


pp.  273-81. 
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Table  6.  Comparison  of  calculated  and  observed  saturation 
magnetostriction  at  15°C. 


Appendix  3. 

A BAND  MODEL  FOR  MAGNETITE 

Throughout  the  preceding  discussion,  we  have  assumed  an  ionic  model 

for  magnetite.  This  approach  provides  an  explanation  of  the  transition  and  of 

the  conductivity.  It  yields  correct  information  about  the  relative  stability  of  the 

normal  and  inverse  spinel  structure  for  the  ferrites,  and  similar  oxides  of 

aluminum  and  chromium/’  and  predicts  on  the  basis  of  a calculation  of  the 

Madelung  energy*  that  the  ordered  structure  of  magnetite  is  slightly  more 

22)  E.  J.  W.  Verwey,  F,  de  Boer  and  J.  H.  var  Santen,  J.  Chem.  Phys. 

16,  1091  (1948). 

* The  actual  value  of  the  calculated  energy  difference  between  the  ordered 
and  the  disordered  structure  is  ca.  1 ev  per  FejO^,  as  compared  to  an 
observed  difference  (from  specific  heat  measurements)  of  0.005  ev,  and 
has  little  significance,  since  it  is  obtained  by  subtracting  two  relatively 
large  numbers.  Also  the  calculated  energy  of  the  ordered  structure  is 
based  on  the  assumption  that  the  ions  occupy  the  same  relative  position# 
in  both  structures,  which  is  not  true;  this  may  account  for  the  discrepancy. 
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stable  then  the  disordered  one.  ' 

Several  attempts  by  Miyahara  2^  to  set  up  a band  model  for  the  ferro- 
magnetic semiconductors  have  not  been  very  successful.  His  basic  assumption 
was  that  a large  exchange  energy  between  electrons  in  the  valence  band  would 
lead  to  ferromagnetism  because  of  holes,  i.e.,  when  electrons  were  removed 
from  the  valence  band  the  remaining  electrons  would  have  a nonzero  spin. 

Any  theory  based  on  this  assumption  must  lead  to  a disappearance  of  the  ferro- 
magnetism at  a sufficiently  low  temperature,  a behavior  thus  far  not  observed 
in  any  material.  Verwey  claims  that  in  the  oxides  of  the  transition  metals,  the 

highest  occupied  levels  are  those  arising  from  metallic  3d  states,  while  the 

25) 

energy  levels  of  the  oxygen  lie  considerably  lower.  ' The  actual  form  of  the 
3d  band  of  the  iron  ions  in  magnetite  would  be  quite  complicated,  since  the 
levels  arising  from  ions  in  octahedral  and  tetrahedral  interstices  will  not  coin- 
cide. The  exchange  energies  responsible  for  the  magnetism  will  cause  compli- 
cated shifts  in  the  levels  depending  on  the  spin  of  the  electron. 

There  is  no  obvious  way  of  estimating  these  various  effects.  Although 
the  whole  picture  is  probably  greatly  oversimplified,  we  can  make  some  ten- 
tative calculations  on  the  following  basis:  the  3d  levels  in  magnetite  are  sup- 
posed to  be  split  into  two  bands,  separated  by  a relatively  narrow  gap.  We 
neglect  the  magnetic  interactions  in  magnetite.  In  distributing  the  electrons 
between  the  two  hands,  we  arrange  that  the  lower  band  shall  be  filled  at  ab- 
solute zero,  and  the  upper  band  empty.  Hence  we  must  have  128  levels/unit 
cell  in  the  valence  band  and  112  levels/unit  cell  in  the  conduction  band.  To 
keep  the  calculations  as  simple  as  possible,  we  have  assumed  that  these  bands 

23)  F.  deBoer,  J.  H.  van  Santen  and  E.  J.  Verwey,  J.  Chem.  Phys.  18, 

1032  (1950).  7 — 

24)  S.  Miyahara,  Phys.  Rev.  55,  105  (1939);  Z.  Physik  113,  247  (1939);  Proc. 
24*56^1942)*  S°C*  Japan"^  49  (1942):  T-  Hirone"an3  S.  Miyahara.  ibid. 

25)  E.  J.  W.  Verwey,  "Semiconducting  Materials,  " Butte rworth,  London,  1951, 
p.  151. 
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are  parabolic.  The  widths  were  obtained  by  requiring  that  /. , (dN/d£)d£ 

equal  the  total  number  of  levels  in  the  band.  At  the  energy  at  which  this 
occurred,  dN/dE  was  assumed  to  drop  discontinuously  to  zero.  This  actually 
had  no  effect  on  the  calculations,  since  at  the  highest  temperature  used,  kT 
was  less  than  one  third  of  the  width  of  the  valence  band.  In  the  calculation  of 
the  number  of  carriers,  we  can  ignore  such  small  effects  and  assume  that  the 
bands  are  parabolic  and  extend  to  infinity. 

There  is  still  nothing  in  this  picture  to  lead  to  the  observed  transition. 

To  account  for  the  transition,  we  suppose  that  each  ferrous  ion  has  a localized 
energy  level,  which  lies  in  the  energy  gap  between  the  two  bands.  Since  the 
energy  gap  itself  is  very  narrow,  these  localized  levels  must  lie  close  to  the 
valence  band.  At  some  temperature  there  will  be  a sufficient  number  of  these 
levels  occupied  so  that  a considerable  overlap  will  occur  between  neighboring 
wave  functions.  When  this  happens  electrons  can  be  transferred  from  one 
localized  level  to  another.  Then  these  levels  will  broaden  slightly  into  a band, 
and  conduction  can  occur  in  this  band. 

To  make  this  idea  quantitative,  we  must  consider  the  distribution  of  the 
ions  in  the  ordered  structure.  Surrounding  a ferrous  ion  is  the  following 
arrangement  of  other  ferrous  ions;  two  at  2.97A;  four  at  5.14A  ; twelve  at 
5.93A;  four  at  6.  63A;  eight  at  7. 85A;  sixat8.39A.  The  radius  of  the  wave 
functions  of  these  localized  levels  will  be  *a0>  where  k is  the  dielectric  con- 
stant and  a is  the  Bohr  radius.  If  we  assume  a dielectric  constant  of  ca.  11, 
o 

the  radius  of  the  wave  function  will  extend  approximately  to  the  twelve  neighbor- 
ing ferrous  ions  at  a distance  of  5.93A,  and  significant  amounts  of  overlap  will 
occur  when  a neighbor  within  a distance  of  5.93A  is  also  excited.  We  expect 

the  transition  to  occur  when  l/l8  of  the  ferrous  ions  are  excited.  Since  there 

20 

are  1.  354  x 1022  ferrous  ions/cc,  the  transition  will  occur  when  7.  52  x 10 
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electrons/cc  have  been  excited  to  local- 


ized levels. 


After  a number  of  trials,  the  band 
scheme  indicated  in  Fig.  20  was  adopted. 
The  effective  masses  for  electrons  and 
boles  are,  respectively,  35  and  55.  The 
widths  of  the  valence  (0.  24  ev)  and  the 
conduction  (0.35  ev)  bands  are  narrow 
compared  to  the  width  of  the  3d  band  in 
the  transition  metals  (ca.  5 ev).  Taking 
the  top  of  the  valence  band  as  the  zero  of 
energy,  we  can  immediately  write 


.35  ev 


112  levels/unit  cell 


8 levels 


.045  ev 


.24  ev 


128  levels 


Fig.  20.  Schematic  diagram  for  the 
energy  levels  of  magnetite. 


p = 1.96xl018T3/2exp(-Ef/kT) 
nc=  9.83z  1017T3/2  exp (Ef  . .12/kt) 
ne=  1.  35xl022(l  + exp (.45  - Ef/kT))_i 


where  p is  the  number  of  holes  in  the  valence  band;  n<;  is  the  number  of  elec- 
trons in  the  conduction  band;  ne  is  the  number  of  electrons  in  the  localized 
levels.  These  equations  for  nc  and  p are  valid  only  if  kT  is  small  compared 
to  E{  and  E{  - 0.12.  Obviously,  these  conditions  will  be  satisfied  only  at  very 
low  temperatures  in  the  present  model.  A convenient  method  of  calculating 
p and  nc  when  kT  is  not  small  was  worked  out  and  used  in  the  calculations. 

This  model  can  account  qualitatively  for  the  observed  axis- switching. 
The  transfer  of  electrons  from  ferrous  to  neighboring  ferric  ions  below  the 
transition  can  occur  through  these  localized  excited  levels.  This  immediately 
leads  u,  to  expect  practically  identical  activation  energies  for  axis-switching 
and  for  the  conductivity  above  the  transition.  The  experimental  values  are 


0.033  and  0.039.  The  band  model  also  accounts  for  the  maximum  observed  in 
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the  conductivity  at  ca.  15°C.  Since  the  number  of  carriers  does  not  increase 
very  much  with  temperature,  the  increased  scattering  at  higher  temperatures 
becomes  sufficient  to  cause  the  conductivity  to  decrease  with  increasing  temper- 
atures. 

In  this  model,  we  have  four  adjustable  parameters:  the  effective  masses 
in  the  valence  and  conduction  bands;  the  width  of  the  gap;  the  location  of  the 
localized  energy  levels  in  the  gap.  There  are  two  conditions  which  we  can  im- 
pose on  this  model.  One  is  that  it  yield  the  correct  transition  temperature  and 
the  other,  that  it  show  the  correct  increase  in  conductivity  at  ca.  550°C.  There 
are  many  combinations  of  the  parameters  which  will  satisfy  these  two  require- 
ments. To  determine  a particular  model,  it  is  necessary  to  include  some  more 
or  less  arbitrary  additional  conditions.  In  the  model  actually  used,  one  of 
these  conditions  was  that  the  Fermi  level  lie  in  the  middle  in  the  middle  of  the 
gap  at  550°C.  Since  the  model  is  so  greatly  over-simplified,  a detailed  com- 
parison with  experimental  data  would  serve  no  purpose.  One  comparison  can 
be  made  with  Hall  effect  measurements  above  the  transition.  Lavine  2^  found 
that  the  effect  in  magnetite  is  so  small  that  it  was  possible  to  determine  an 
upper  limit  only  for  the  Hall  coefficient,  i.e.,  a lower  limit  for  the  number  of 
carriers.  Assuming  that  only  the  electrons  contribute  to  the  conductivity, 

Lavine  has  calculated  that  the  density  must  be  greater  than  102*/cc.  This  is  in 
agreement  with  the  model,  which  predicts  values  of  n^  from  1.  54  x 102^  at  -99° 
to  3.66  x 1021  at  17°C. 

This  model  is  satisfactory  as  a first  and  very  crude  approximation  to 
the  band  structure  of  magnetite.  The  present  model  is  basically  similar  to 
the  ionic  order -disorder  picture. 

26)  L.  M.  Lavine,  Prog.  Rep.  26,  Cruft  Laboratory,  Harvard  University, 

Jan.  , 1953,  p.  32. 


